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In this work we reconsider the problem discussed by Caldeira and Leggett (CL), A. O. Caldeira 
and A. J. Leggett, Ann. Phys., 149, 374 (1983), and generalize their results for a quantum oscillator 
^ , s coupled bilinearly to a reservoir R with dense discrete spectrum of harmonic oscillators with 

close to uj 3 frequencies. We show that for such systems dynamic evolution of any state of the CL 
Hamiltonian consists of recurrence cycles with partial revivals of the initial state amplitude. This 
revival or Loschmidt echo appears in each cycle due to reverse R to s transitions (not existing in 

■ the limit of a macroscopic bath). Width and number of the Loschmidt echo components increase 
with the recurrence cycle number fc, eventually leading to irregular, stochastic-like time evolution. 
Standard for continuous spectrum thermal bath CL dynamics dissipative dynamics (exponential 

■ decay) takes place only within the initial cycle k = 0. In terms of the effective CL Hamiltonian, 
where the reservoir degrees of freedom are integrated out, our results mean oscillatory (time periodic) 
behavior of the effective mass. Only in the initial cycle, we found the known for CL Hamiltonian 

I ■ exponential in time increase of the effective mass. We do believe that the basic ideas inspiring our 

work can be applied to a large variety of interesting nano-systems and molecular complexes. 

PACS numbers: 03.65, 82.20.B, 05.45.-a, 72.10.-d 
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Introduction. In this work we are dealing with dynamics of nano-particles or relatively large molecular complexes. 
In this case the system has a large (of the order of 10 2 — 10 4 ) number of internal vibration modes with a dense discrete 
spectrum. Characteristic mean interlevel spacing Ao ~ 1 — 100 cm -1 , in high-frequency region (1000 — 2000 cm" 1 ) 
corresponds to Poincare recurrence periods in the range 0.1 — IQps routinely probed by modern optical methods. 
Such dynamics of the initial excited state coupled to the reservoir of vibration states is in the heart of the realm of 
radiationless transitions [l| , Q , , 0] ■ These experiments show that in many cases characteristic life time of lowest 
■ excited states, is on the order of recurrence cycle period. Therefore the both phenomena (namely, decay of the state s 
and exchange with R states) has to be treated together. Although all robust dynamic characteristics of radiationless 
transitions are described by the CL Hamiltonian [f|, @ which reads (in the self-evident notations with p s and q s 
\& \ momentum and coordinate of the harmonic oscillator coupled bilinearly (with coupling C n ) to the discrete states of 
the reservoir, also formed by harmonic oscillators p n and q n ) 

H=\(P 2 S + uUl) + \ Y^ivl + "111) (1) 

n n 

the Hamiltonian (JTJ has been investigated with all details only in the limit of continuous spectrum thermal bath. In 
this limit the stored in the excited state energy is exponentially in time dissipated via Ohmic friction. Therefore many 
relevant for nano-systems and relatively large molecular complexes dynamic aspects of revivals of the initial states 
are merely disregarded. Surprisingly for us, scanning the literature we did not find any paper treating theoretically 
discrete spectrum CL model. 

Our model is applicable to a standard setup used for radiationless transitions, and includes the following main 
ingredients. Preparation of the initial excited state s, coupled strongly (matrix elements C n ) to other high frequency 
vibration states {n}. These states form reservoir R with a dense discrete spectrum of levels distributed around the 
excited initial state, i.e., resonance region width is smaller than characteristic frequencies 

\u - lu s \ < uj s , u„ . (2) 
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In the continuous spectrum reservoir one can integrate out the reservoir degrees of freedom and then end up with the 
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following effective Hamiltonian Fl%tt bearing dissipative physics Q 

H c eff = eM-Vt)\Ps + cM+rt)\Ql , (3) 

where entering effective Hamiltonian canonical momentum and coordinate P s and Q s (satisfying the standard com- 
mutation relations), and "friction" coefficient 77 characterizes the coupling to the reservoir, and it determines the time 
dependent effective mass. Eq. ([3]) is a Hamiltonian of damped oscillator under external force determined by initial 
values of the reservoir momenta and coordinates. Irreversible transition to the reservoir states yields to exponential 
decay of the initial state, and that is all what one can expect from full CL (TTJ) (or its reduced version Hamil- 
tonian for the continuous thermal bath. The aim of our work is to see how this conclusion above (and the effective 
Hamiltonian ([5])) has to be modified for systems where the vibrational state reservoir is characterized by discrete 
spectrum. 

Secular equation for CL Hamiltonian with discrete spectrum reservoir. The equations of motion for 
the Hamiltonian of bilinearly coupled oscillators (TTJ 

q s + us 2 = - ^2 C nq n ; q n + u\q n = -C n q s , (4) 

n 

supplemented by the initial conditions 

q s (0) = q s a , q s {0) = q s a , q n {0) = q n o , <M) (5) 
can be easily solved by Laplace transformation Q , [§| 



<7s,{n}(w) = / q s ,{n}(t) exp(ut)dt. (6) 
Jo 

Then the secular equation to find eigen-frequencies 

FH^ 2 -^ 2 ^^^. (7) 

n n 

In the resonance conditions ©, Eq. ([7]) reads as 

F{e) = F (e) + eF x {e) , (8) 

where e = u) — uj s and 

F (e) = 2. s L - £ JL) ; Fl (e) = e + £ ( 1 + £■) ' , (9) 



with the renormalized coupling constants 



Zwanzig approximation. To simplify the calculations we start with Zwanzig model |8J assuming equidistant in 
the resonance region the reservoir spectrum uj n = lo s + n and B„ = B. Then the secular equation @ acquires the 
compact form 

F (e) = e-77C0t(7re),, (11) 

where the friction coefficient 77 = ttB 2 is the full rate of the transitions from the s oscillator into all R states In terms 
of the harmonic oscillator coupled to a discrete spectrum reservoir, its behavior is dramatically different from that 
for the continuous spectrum thermal bath. In the latter case the friction coefficient ?y is a time independent constant. 
Therefore the coupling yields to irreversible decay (dissipation of the initially excited state). It is not the case for the 
discrete spectrum. The effective "friction coefficient" is time dependent (e.g., it equals zero in a certain time within 
each of the recurrence cycle), and it leads to Loschmidt echo phenomenon and non-trivial quantum dynamics. 
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CL effective Hamiltonian for discrete spectrum reservoir. Having in hands the CL Hamiltonian and the 
equations of motion we can express the amplitude of the s oscillator as a series over the initial values with time 
dependent coefficients 

q s (t) = oi(t)g,(0) + oa(t)g.(0) + X)(°ni(*)*»( ) + a n2(t)q n {0)) , (12) 

n 

and similarly for the Laplace transforms 

q s (u) = ai(u)q s (0) + a 2 {u)q s (0) + ^(a nl (u)q n (0) + a n2 (u)q n (0)) . (13) 

n 

This expansion (fT2]l is determined by the poles of the function l/i*b(e) which are the roots of the secular equation. 
Namely the time dependent coefficients are 

ai(i) - (cos(w.t) + -5- smfat)) W(t) ; a 2 (i) = SmM W(t) . (14) 

Slow (in a scale w^ 1 ) variations of the amplitude of s oscillator is described by the function W, which is represented 
in a Fourier series over the Hamiltonian eigen frequencies 

ur(+\ ST 27?cos(e»t) 



n>0 



The series is calculated by Poisson summation formula and then it is reduced to a sum over partial recurrence cycle 
amplitudes H, 03, El 

[t/2n] 

W(t) = J2 i W ^ = -? exp(-r fc /2)0(r fc )4_ 1 (r fc ) , (16) 

fe=0 

where, [x] means the integer part of x, 0(x) is the Heaviside step-function, Lj._ 1 is adjoint Laguerre polynomial [l2j, 
and the local recurrence cycle "time" is defined as = r)(t — 2kn). These expressions (IT21 - (|TE)) manifest that time 
evolution of the excited initial state coupled to the discrete spectrum reservoir can be reduced to one dimensional 
damped oscillator with time dependent friction. Only for the initial k = cycle, when from (fl"6]l Wo(t) = exp(— r]t/2) 
we get simple exponential decay as it is the case for CL dynamics with the continuous spectrum thermal bath. From 
(|12p. (|16[) we conclude that time evolution of q s and p s is determined by the functions a\(t) and a-zit) which depend 
on the oscillator frequency uj s and friction coefficient r\ but independent of time evolution of the reservoir states in a 
phase space {q n , p n } orthogonal to the phase space q s , p s . Then it is convenient to make a projection Q 



and construct the monodromy matrix 



,(t) = oi(t)«.(0) + o 2 (t)«.(0), (17) 
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The determinant of the monodromy matrix (|18p equals to VF 2 (t) and then following the known procedure [L 
we end up with the generalized for discrete reservoir spectrum effective CL Hamiltonian 



H d eff = W\t)\p 2 8 +W-\t)^Ql. (i<>) 



The derivation of this effective Hamiltonian Hf^ is the main result of our paper and it is ready for further inspection 
and various applications. Unlike the standard CL Hamiltonian ([2]) where with time going on the oscillator s irreversibly 
loses its energy and therefore the amplitude of the oscillations decays exponentially, for our Hamiltonian (fT9|) this 
picture holds only for the initial k — cycle. Then the oscillator s not only transfers its energy to the discrete 
spectrum R oscillators but as well receives the energy from the reservoir. We show time evolution of the oscillator s 
amplitude in the fig. 1. In the cycle k = the initial s-oscillator state decays exponentially (the upper panel in the 
fig. 1). This decay is accompanied by the monotonous growth of the R oscillators amplitudes (second and third from 
the top panels in the fig. 1). In the cycle k = 1 the s oscillator amplitude increases because of the reverse R s 
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transitions. Meanwhile the states of the reservoir are depopulated-populated almost synchronously and it yields to the 
intense one-component Loschmidt echo signal. In the next k > 1 cycles, the synchronisation of the reverse transition is 
gradually destroyed. Then the Loschmidt echo signals become multi-component ones. Eventually, when the width of 
the echo signal becomes comparable with the cycle period, one observes mixing phenomena and crossover to irregular, 
stochastic like dynamics [lH, [l6j . 

Evolution of quantum states of the s oscillator. To get the standard quantum mechanical coordinate - 
momentum commutation relation we have to dehne the momentum P s conjugated to Q s (|17p as an operator 

P s = W(t)Q s . (20) 

The same function W(t) determines non-zero matrix-elements 

{Q% v -i = V -W 2 {t). (21) 

Therefore if the s oscillator initial (t = 0) state has the quantum number A, the wave function ^A(Q s ,t) a t the time 
t is uniquely determined by the eigen function $a(Qs) of the initial state 

* A (Qs,t)=A s (A,t)$ A (Q s ), (22) 

where the wave function amplitude A s (A,t) = (W(t)) A . We conclude from the ([22]) that the ground state wave 
function is time independent (as it should be), whereas the wave functions for the excited states change in time 
manifesting recurrence cycle behavior. For the initial cycle k = 

A S (A, t) = exp(-A??t/2) . (23) 

Because in the harmonic approximation only transitions where the quantum numbers are changed by ±1 are not for- 
bidden, the global (s oscillator and R oscillators) wave function is the superposition (with time dependent coefficients) 
of the harmonic oscillator wave functions with the total quantum number A, i.e., 

* A (t)= A(v s ;{v n };t)$ Us (q s ) ( 24 ) 

v s ;{i>n} n 

where mentioned above selection rule requires that v s + ^2 n v n = A. We show in the Fig. 2, the time evolution of 
the A = 4 initial state of the s oscillator. For the excited state, the superposition includes the wave function of 

the ground state (A s {0„}) ; ((A — l) s 1„ {0 n '=£n}) The coefficients A entering the ([24| are determined by 

the product of the corresponding powers of the function W(t). The echo signal decreases with A, and as a result, 
the decay of the highly excited initial state becomes irreversible (the upper panel in the fig. 2). Then less excited 
states also decay, and in the next cycles, the excitation energy is redistributed (as one can expect approaching to the 
equilibrium state) between low energy excited states of the R oscillators (the lowest panel in the fig. 2). 

Conclusion. To summarize, in this paper we investigated quantum dynamics for a quantum oscillator s coupled 
bilinearly to a reservoir R with dense discrete spectrum of harmonic oscillators with close to uj s frequencies. We 
show that for such systems dynamic evolution of any state of the CL Hamiltonian is reduced to recurrence cycles 
(oscillations) due to reverse (from the discrete spectrum reservoir R to s) transitions not existing in the limit of 
a macroscopic thermal bath. These reverse transitions lead to a partial restoration of the s oscillator amplitudes 
(Loschmidt echo). Width and number of the Loschmidt echo components increase with the recurrence cycle number 
k, eventually leading to irregular, stochastic-like time evolution. As a possible physical realization and evidence for 
the discussed above phenomena we note recent experiments on ultra-fast (in ps scale) energy exc hang e between high 



frequency vibrations 17| and observations of almost ballistic energy transport in molecular chains |18[ . We anticipate 
also that considered above behavior of the CL Hamiltonian is relevant to describe experimental data for radiationless 
transitions and ultrafast photo-chemical reactions. 
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Figure captions 

Fig. 1 

Time evolution of the oscillator s state with A = 1, and for the R oscillators n = 1 , 3 , 4 , 7 which were initially in 
the ground state. The friction coefficient r\ — 2.2. 

Fig. 2 

Time evolution of the oscillator s state with A = 4, and occupations of the following states s and R oscillators: 
(3 S 1„) , (2 S 1„ \ n < , (l s 1„ 1„/ l n «) , (1„ 1„/ 1„» !„/»)). The friction coefficient r\ = 3.1. 
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